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Abstract 

We present a likelihood analysis of the parameter space of the constrained minimal super- 
symmetric extension of the Standard Model (CMSSM), in which the input scalar masses tuq 
and fermion masses mi/2 are each assumed to be universal. We include the full experimental 
likelihood function from the LEP Higgs search as well as the likelihood from a global precision 
electroweak fit. We also include the likelihoods for 6 — 57 decay and (optionally) (7^ — 2. For 
each of these inputs, both the experimental and theoretical errors are treated. We include 
the systematic errors stemming from the uncertainties in rrit and ruh, which are important 
for delineating the allowed CMSSM parameter space as well as calculating the relic density 
of supersymmetric particles. We assume that these dominate the cold dark matter density, 
with a density in the range favoured by WMAP. We display the global likelihood function 
along cuts in the (mi/2,mo) planes for tan/5 = 10 and both signs of fi, tan/3 = 35,/i < 
and tan/5 = 50, yU > 0, which illustrate the relevance oi — 2 and the uncertainty in m^. 
We also display likelihood contours in the (^«i/2,'n^o) planes for these values of tan/5. The 
likelihood function is generally larger for /i > than for /i < 0, and smaller in the focus-point 
region than in the bulk and coannihilation regions, but none of these possibilities can yet be 
excluded. 



CERN-TH/2003-262 
October 2003 



1 Introduction 



Supersymmetry remains one of the best-motivated frameworks for possible physics beyond 
the Standard Model, and many analyses have been published of the parameter space of the 
minimal supersymmetric extension of the Standard Model (MSSM) . It is often assumed that 
the soft supersymmetry-breaking mass terms mi/2,mo are universal at an input GUT scale, 
a restriction referred to as the constrained MSSM (CMSSM). In addition to experimental 
constraints from sparticle and Higgs searches at LEP [1] , the measured rate for 6 — > 57 [2] and 
the value of — 2 [3] ^, the CMSSM parameter space is also restricted by the cosmological 
density of non-baryonic cold dark matter, flcDM [5, 6, 7, 8]. It is also often assumed that 
most of flcDM is provided by the lightest supersymmetric particle (LSP), which we presume 
to be the lightest neutralino x- 

The importance of cold dark matter has recently been supported by the WMAP Collab- 
oration [9, 10], which has established a strong upper limit on hot dark matter in the form of 
neutrinos. Moreover, the WMAP Collaboration also reports the observation of early reion- 
ization when 2 ~ 20 [10], which disfavours models with warm dark matter. Furthermore, 
the WMAP data greatly restrict the possible range for the density of cold dark matter: 
^CDuh^ — 0.1126lo;oo9i (one-cr errors). Several recent papers have combined this infor- 
mation with experimental constraints on the CMSSM parameter space [11, 12, 13, 14, 15], 
assuming that LSPs dominate flcDM- 

The optimal way to combine these various constraints is via a likelihood analysis, as has 
been done by some authors both before [16] and after [13] the WMAP data was released. 
When performing such an analysis, in addition to the formal experimental errors, it is also 
essential to take into account theoretical errors, which introduce systematic uncertainties 
that are frequently non-negligible. The main aim of this paper is to present a new likelihood 
analysis which includes a careful treatment of these errors. 

The precision of the WMAP constraint on Qcdm selects narrow strips in the CMSSM 

parameter space, even in the former 'bulk' region at low mi/2 and niQ. This narrowing is 

even more apparent in the coannihilation 'tail' of parameter space extending to larger nii/2, 

in the 'funnels' due to rapid annihilations through the A and H poles that appear at large 

tan/3, and in the focus-point region at large rriQ, close to the boundary of the area where 

elcctrowcak symmetry breaking remains possible. The experimental and theoretical errors 

are crucial for estimating the widths of these narrow strips, and also for calculating the 

^In view of the chequered history of this constraint, we present results obtained neglecting 5^ — 2, as well 
as results using the latest re-evaluation of the Standard Model contribution [4] . 
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likelihood function along cuts across them, as well as for the global hkelihood contours we 
present in the {mi/2,mo) planes for different choices of tan/3 and the sign of /i. 

In the 'bulk' and coannihilation regions, we find that the theoretical uncertainties are 
relatively small, though they could become dominant if the experimental error in ilcDuh^ 
is reduced below 5% some time in the future. However, theoretical uncertainties in the 
calculation of ruh do have an effect on the lower end of the 'bulk' region, and these are 
sensitive to the experimental and theoretical uncertainties in rrit and (at large tan/3) also 
nib. The theoretical errors due to the current uncertainties in and rrit are dominant in 
the 'funnel' and 'focus-point' regions, respectively. These sensitivities may explain some of 
the discrepancies between the results of different codes for calculating the supersymmetric 
relic density, which are particularly apparent in these regions. These sensitivities imply that 
results depend on the treatment of higher-order effects, for which there are not always unique 
prescriptions. 

With our treatment of these uncertainties, we find that the half-plane with // > is 
generally favoured over that with // < 0, and that, within each half- plane, the coannihilation 
region of the CMSSM parameter space is generally favoured over the focus-point region ^, 
but these preferences are not strong. 

The rest of the paper is organized as follows. In section 2, we discuss the treatment of 
the various constraints employed to define the global likelihood function. In section 3, we 
present the profile of the global likelihood function along cuts in the (mi/2,mo) plane for 
different choices of tan /3 and the sign of In section 4, we present iso-hkelihood contours 
at certain CLs, obtained by integrating the likelihood function. Finally, in section 5, we 
summarize our findings and suggest directions for future analyses of this type. 

2 Constraints on the CMSSM Parameter Space 
2.1 Particle Searches 

We first discuss the implementation of the accelerator constraints on CMSSM particle masses. 
Previous studies have shown that the LEP hmits on the masses of sparticles such as the 
selectron and chargino constrain the CMSSM parameter space much less than the LEP 
Higgs hmit and b ^ s'y (see, e.g., [7, 17]). As we have discussed previously, in the CMSSM 
parameter regions of interest, the LEP Higgs constraint reduces essentially to that on the 
Standard Model Higgs boson [17]. This is often implemented as the 95% confidence-level 
^Our conclusions differ in this respect from those of [13]. 
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lower limit ruh > 114.4 GeV [1]. However, here we use the full hkehhood function for the 
LEP Higgs search, as released by the LEP Higgs Working Group. This includes the small 
enhancement in the likehhood just beyond the formal hmit due to the LEP Higgs signal 
reported late in 2000. This was re-evaluated most recently in [1], and cannot be regarded 
as significant evidence for a light Higgs boson. We have also taken into account the indirect 
information on rrih provided by a global fit to the precision electroweak data. The likelihood 
function from this indirect source does not vary rapidly over the range of Higgs masses found 
in the CMSSM, but we include this contribution with the aim of completeness. 

The interpretation of the combined Higgs likelihood, Cexp, in the (mi/2, mo) plane depends 
on uncertainties in the theoretical calculation of m^. These include the experimental error 
in rrit and (particularly at large tan /3) nib, ^-nd theoretical uncertainties associated with 
higher-order corrections to nih. Our default assumptions are that = 175 ± 5 GeV for the 
pole mass, and nib = 4.25±0.25 GeV for the running MS mass evaluated at nib itself. The 
theoretical uncertainty in nih, (7th i is dominated by the experimental uncertainties in nit^bi 
which are treated as uncorrelated Gaussian errors: 



The Higgs mass is calculated using the latest version of FeynHiggs [18]. Typically, we find 
that (dmh/dnit) ~ 0.5, so that ath is roughly 2-3 GeV. Subdominant two-loop contributions 
as well as higher-order corrections have been shown to contribute much less [19]. 

The combined experimental likelihood, Cg^p, from direct searches at LEP 2 and a global 
electroweak fit is then convolved with a theoretical likelihood (taken as a Gaussian) with 
uncertainty given by ath from (1) above. Thus, we define the total Higgs likelihood function, 
£h, as 



where jV" is a factor that normalizes the experimental likelihood distribution. 
2.2 b ^ s'j Decay 

The branching ratio for the rare decays 6 — > 57 has been measured by the CLEG, BELLE 
and BaBar collaborations [2], and we take as the combined value B{b — > S7) = (3.54 ± 0.41 ± 
0.26) X 10~^. The theoretical prediction of 6 — > 57 [20, 21] contains uncertainties which stem 
from the uncertainties in nib, o^s, the measurement of the semileptonic branching ratio of the 
B meson as well as the effect of the scale dependence. In particular, the scale dependence 
of the theoretical prediction arises from the dependence on three scales: the scale where the 




(1) 
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QCD corrections to the semileptonic decay are calculated and the high and low energy scales, 
relevant to 6 ^ 57 decay. These sources of uncertainty can be combined to determine a total 
theoretical uncertainty. Finally, the experimental measurement is converted into a Gaussian 
hkelihood and convolved with a theoretical hkelihood to determine the total hkelihood Cijgg 
containing both experimental and theoretical uncertainties [20] ^. 

2.3 Measurement oi — 2 

The interpretation of the BNL measurement of = (7^ — 2 [3] is not yet settled. Two 
updated Standard Model predictions for have recently been calculated [4]. One is based 
on e~^e~ hadrons data, incorporating the recent re-evaluation of radiative cross sections 
by the CMD-2 group: 

Uf, = (11, 659, 180.9 ± 7.2 ± 3.5 ± 0.4) x 10~^°, (3) 

and the second estimate is based on r decay data: 

= (11, 659, 195.6 ± 5.8 ± 3.5 ± 0.4) x 10"^°, (4) 

where, in each case, the first error is due to uncertainties in the hadronic vacuum polarization, 
the second is due to hght-by-hght scattering and the third combines higher-order QED and 
electroweak uncertainties. Comparing these estimates with the experimental value [3], one 
finds discrepancies 

Aa^ = (22.1 ± 7.2 ± 3.5 ± 8.0) x 10-^°(1.9 a) (5) 

and 

Aa^ = (7.4 ± 5.8 ± 3.5 ± 8.0) x 10-^°(0.7 a), (6) 

for the e'^e~ and r estimates, respectively, where the second error is from the light-by-light 
scattering contribution and the last is the experimental error from the BNL measurement. 

Based on the e~^e~ estimate, one would tempted to think there is some hint for new 
physics beyond the Standard Model. However, the r estimate does not confirm this optimistic 
picture. Awaiting clarification of the discrepancy between the e'^e~ and r data, we calculate 
the likelihood function for the CMSSM under two hypotheses: 

• neglecting any information from gf^ — 2, which may be unduly pessimistic, and 

^Further details of our treatment of experimental and theoretical errors, as applied to the CMSSM, can 
be found in [5]. 
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• taking the e'^e estimate (5) at face value, which may be unduly optimistic. 



When including the likelihood for the muon anomalous magnetic moment, a^, we calcu- 
late it combining the experimental and the theoretical uncertainties as follows: 

where cr^ = al^p + cr^^, with Uexp taken from (5) and crf^ from (1), replacing rrih by a^. 

As is well known, the discrepancy (5) would place significant constraints on the CMSSM 
parameter space, favouring > 0, though we do consider both signs of In fact, we 
find that > is favoured somewhat, even with the 'pessimistic' version (6) of the — 2 
constraint. 

2.4 Density of Cold Dark Matter 

As already mentioned, we identify the relic density of LSPs with QcDhih^- In addition to the 
CMSSM parameters, the calculation of VtcDuh^ involves some parameters of the Standard 
Model that are poorly known, such as and m^. The default values and uncertainties we 
assume for these parameters have been mentioned above. Here we stress that both these 
parameters should be allowed to run with the effective scale Q at which they contribute 
to the calculation of the relic density, which is typically Q ~ 2m^. This effect is particu- 
larly important when treating the rapid-annihilation channels due to XX A, H XX 
annihilations, but is non-negligible also in other parts of the CMSSM parameter space. 

Specifically, the location of the rapid-annihilation funnel due to A, H Higgs-boson ex- 
change, which appears in the region where uia — 2m^, depends significantly on the deter- 
mination of rriA [6]. For this determination, the input value of the running MS mass of rrib 
is a crucial parameter, and the appearance of the funnels depends noticeably on [5, 22]. 
On the other hand, the exact location of the focus-point region [23] (also known as the hy- 
perbolic branch of radiative symmetry breaking [24]) depends sensitively on rrit [25, 22, 7], 
which dictates the scale of radiative electroweak symmetry breaking [26]. 

In calculating the likelihood of the CDM density, we follow a similar procedure as for the 
anomalous magnetic moment of the muon in (1, 7), again taking into account the contribution 
the uncertainties in rrit^b- In this case, we take the experimental uncertainty from WMAP 
[9, 10] and the theoretical uncertainty from (1), replacing rrih by VL^h"^. We will see that the 
theoretical uncertainty plays a very significant role in our analysis. 
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2.5 The Total Likelihood 



The total likelihood function is computed by combining all the components described above: 

C'tot — ChX Cbsj X C(^^h'^ X (8) 

In what follows, we consider the CMSSM parameter space at fixed values of tan (3 — 10, 35, 
and 50 with Aq — 0. For tan/3 = 10 and 35, we compute the hkelihood function for both 
signs of but not for tan/? = 50, since in the case the choice // < does not provide a 
solution of the RGEs with radiative clcctroweak symmetry breaking. 

The likelihood function in the CMSSM can be considered as a function of two variables, 
J^toti'm'i/2,'>TiQ), where mi/2 and mo are the unified GUT-scale gaugino and scalar masses 
respectively. When plotting confidence levels as iso-likelihood contours in the corresponding 
(mi/2,mo) planes, we normalize likelihood function by setting the volume integral 

j Ctotdmodmi/2 = 1 (9) 

for each value of tan (3, combining where appropriate both signs of /i. We also compare the 
integrals of the likelihood function over the coannihilation and focus-point regions, and for 
different values of tan p. 

For most of the results presented below, we perform the analysis over the range mi/2 = 
100 GeV up to 1000 GeV for tan /3 = 10 and up to 2 TeV for tan /3 = 35 and 50. The upper 
limit on mo is taken to be the limit where solutions for radiative electroweak symmetry 
breaking are possible and the range includes the focus-point region at large mo- We discuss 
below the sensitivity of our results to the choice of the upper limit on mi/2- 

3 Widths of Allowed Strips in the CMSSM Parameter 
Space 

We begin by first presenting the global likehhood function along cuts through the (mi/2, mo) 
plane, for different choices of tan/3, the sign of /i and mi/2. These exhibit the relative 
importance of experimental errors and other uncertainties, as well as the potential impact 
of the Qfj, — 2 measurement. 

We first display in Fig. 1 the hkelihood along slices through the CMSSM parameter 
space for tan/3 = 10,^0 = 0, // > 0, and mi/2 = 300 and 800 GeV in the left and right 
panels, respectively, plotting the likelihood as a function of mo in the neighborhood of the 
coannihilation region [27]. The sohd red curves show the total likehhood function calculated 
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Figure 1: The likelihood function along slices in rriQ through the CMSSM parameter space 
for tan/5 = IOjAq = 0,/x > and mi/2 = 300,800 GeV in the left and right panels, 
respectively. The solid red curves show the total likelihood function and the green dashed 
curve is the likelihood function with Arrit = Ami, = 0. Both analyses include the — 2 
likelihood calculated using e~^e~ data. The horizontal lines show the 68% confidence level of 
the likelihood function for each case. 

including the uncertainties which stem from the experimental errors in mt and m;,. The 
green dashed curves show the likelihood calculated without these uncertainties, i.e., we set 
Am-f = Am^ = 0. We see that these errors have significant effects on the likelihood function. 
In each panel, the horizontal lines correspond to the 68% confidence level of the respective 
likelihood function. The likelihood functions shown here include Ca^ calculated using e^e~ 
data. For these values of mi/2 and with /i > 0, the constraint from g^ — 2 is not very 
significant. For reference, we present in Table 1 and 2 the values of the likelihood functions 
corresponding to the 68%, 90%, and 95% CLs for each choice of tan/3 and Amj. 

When /i < 0, the g^ — 2 information plays a more important role, as exemplified in 
Fig. 2, where we show the likelihood in the coannihilation region for mi/2 = 800 GeV. For 
mi 12 = 300 GeV, the likelihood is severely suppressed (see the discussion below) and we do 
not show it here. 

We now discuss the components of the likelihood function which affect the relative heights 
along the peaks shown in Fig. 1. In the case mi/2 = 300 GeV, the likelihood increases when 
the errors in mt and mb are included, due to two dominant effects. 1) The total integrated 
likelihood is decreased when the errors are turned on (by a factor of ~ 2 when (yf^ — 2 is 
included and by a factor of ~ 3 when it is omitted, for tan/3 = 10), so the normalization 
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Table 1: Likelihood values (xlO^), including — 2, for the 68%, 90%, and 95% CLs for 
different choices of tan f3 and the uncertainty in rrit. 



tan/3 


CL 


Amt = 5 GeV 


Amt = 1 GeV 


Amt = 0.5 GeV 


Amt = GeV 




68% 


0.14 


0.13 


0.087 


0.046 


10 


90% 


3.0 xlO-3 


1.4 xlO-'^ 


2.5 xlO"^ 


0.021 




95% 


2.9 xlO-5 


6.2 xlO-5 


1.1 xlO-4 


0.011 




68% 


2.2 xlO-* 


2.0 xlO'^ 


2.7 xlO-^ 


1.1x10^ 


35 


90% 


2.8 xlO"^ 


5.0 xlO^^ 


7.5 xlO"^ 


1.8x10^^ 




95% 


1.1 xlO"^ 


2.7 xlO^^ 


3.9 xlO^^ 


6.8x10^^ 




68% 


5.3 xlO-4 


5.7 xlO-^ 


5.4 xlO-^ 


7.0 xl0^4 


50 


90% 


8.7 xlO-5 


7.7 xlO-5 


1.0 xlO-^ 


1.9 xlO-^ 




95% 


2.3 xlO-^ 


3.2 xlO-^ 


4.6 xlO-^ 


8.2 xlO-^ 



Table 2: Likelihood values (xlO^), excluding g^ — 2, for the 68%, 90%, and 95% CLs for 
different choices o/tan/3 and the uncertainty in mt- 



tan/3 


CL 


Amt = 5 GeV 


Amt = 1 GeV 


Amt = 0.5 GeV 


Amt = GeV 




68% 


0.059 


1.6 xlO-^ 


9.6 xlO-^ 


0.052 


10 


90% 


5.6 xlO-5 


1.3 xlO"^ 


2.4 xlO"^ 


0.024 




95% 


4.0 xlO--' 


7.9 xlO-5 


1.3 xlO-'^ 


0.011 




68% 


2.4 xlO-^ 


1.9 xlO"^ 


2.4 xlO-'^ 


7.8x10-^ 


35 


90% 


2.3 xlO-5 


5.0 xlO-5 


7.7 xlO-5 


1.6x10"^ 




95% 


1.2 xlO-^ 


2.8 xlO-^ 


4.0 xlO-^ 


7.5x10-5 




68% 


3.3 xlO-* 


3.5 xlO-^ 


4.3 xlO-* 


6.4 xlO-* 


50 


90% 


4.2 xlO"^ 


6.4 xlO-5 


1.0 xlO"^ 


1.8 xlO"^ 




95% 


2.0 xlO"^ 


3.4 xlO"^ 


5.0 xlO"^ 


8.5 xlO-5 
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Figure 2: .4s %n Fig. 1 for tan/3 = 10, = 0, yU < and mi/2 = 800 GeV. The g^, - 2 
constraint is included (excluded) in the left (right) panels. In the right panel the 68% CLs 
for both cases are incidentally closed to each other. 

constant, A/", becomes larger, and 2) since mi/2 = 300 GeV corresponds to the lower limit 
on mi/2 due to the experimental bound on the Higgs mass, the Higgs contribution to the 
likelihood increases when the uncertainties in the heavy quark masses are included. When 
= 800 GeV, it is primarily the normalization effect which results in an overall increase. 
The Higgs mass contribution at this value of mi/2 is essentially Ch^^^ = 1. We remind the 
reader that the value of the likelihood itself has no meaning. Only the relative likelihoods 
(for a given normalization) carry any statistical information, which is conveyed here partially 
by the comparison to the respective 68% CL likelihood values. 

In Fig. 3, we extend the previous slices through the CMSSM parameter space to the 
focus-point region at large niQ. The solid (red) curve corresponds to the same likelihood 
function shown by the solid (red) curve in Fig. 1, and the peak at low mo is due to the 
coannihilation region. The peak at mo — 2500(4500) GeV for mi/2 = 300(800) GeV is due 
to the focus-point region ^. The g^ — 'i constraint is not taken into account in the upper two 
figures of this panel. Also shown in Fig. 3 are the 68%, 90%, and 95% CL lines, corresponding 
to the iso-likelihood values of the fully integrated likelihood function corresponding to the 
solid (red) curve. As one can see, one of the effects of the g^ — 2 constraint (even at its 
recently reduced significance) is a suppression of the hkelihood function in the focus-point 
region. 

"^We should, in addition, point out that different locations for the focus-point region are found in different 
theoretical codes, pointing to further systematic errors that are currently not quantifiable. 
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m„(GeV) mo(GeV) 



Figure 3: As in Fig. 1, hut for slices at fixed mi/2 that include also the focus-point region at 
large mo. The red (solid) curves are calculated using the current errors in rrit and mi,, the 
green dashed curve with no error in mt, the violet dotted lines with Amt = 0.5 GeV, and the 
blue dashed-dotted lines with Amt = 1 GeV. In the upper two figures, the Qfi — 2 constraint 
has not been applied. 
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m„ (GeV) 



Figure 4: The value of Q^h"^ (solid) and dQh'^/dmt (dashed) as functions of mo /or tan /3 = 
10,y4o = 0,yU > and mi/2 = 300 GeV, corresponding to the slice shown in Fig. 3c. 

The focus-point peak is suppressed relative to the coannihilation peak at low mo because 
of the theoretical sensitivity to the experimental uncertainty in the top mass. We recall that 
the likelihood function is proportional to cr^^, and that a which scales with d{Q^h'^)/dmt, 
is very large at large mo [22]. This sensitivity is shown in Fig. 4, which plots both Q^h"^ and 
d{fl^h'^) / drrit for the cut corresponding to Fig. 3c. Notice that, for the two values of mo 
with ^x^'^ ~ ^-^^ corresponding to the coannihilation and focus-point regions, the error due 
to the uncertainty in m^ is far greater in the focus-point region than in the coannihilation 
region. Thus, even though the exponential in Cn^h'^ is of order unity near the focus-point 
region when fi^/i^ ^0.1, the prefactor is very small due the large uncertainty in the top 
mass. This accounts for the factor of > 1000 suppression seen in Fig. 3 when comparing the 
two peaks of the solid red curves. 

We note also that there is another broad, low- lying peak at intermediate values of mo. 
This is due to a combination of the effects of cr in the prefactor and the exponential. We 
expect a bump to occur when the Gaussian exponential is of order unity, i.e., Q^h'^ ~ 
\/2Amtdilxh'^ /drrit. From the solid curve in Fig. 4, we see that il^h'^ ~ 10 at large mo for 
our nominal value rrit = 175 GeV, but it varies significantly as one samples the favoured 
range of rrit within its present uncertainty. The competition between the exponential and 
the prefactor would require a large theoretical uncertainty in dQ^h'^/drrit ~ 2 for 

Arrit = 5 GeV. From the dashed curve in Fig. 4, we see that this occurs when mo ~ 1000 
GeV, which is the position of the broad secondary peak in Fig. 3a. At higher mo, cr continues 
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to grow, and the pref actor suppresses the hkehhood function until Jl-^/i^ drops to ~ 0.1 in 
the focus-point region. 

As is clear from the above discussion, the impact of the present experimental error in rrit 
is particularly important in this region. This point is further demonstrated by the differences 
between the curves in each panel, where we decrease ad hoc the experimental uncertainty in 
rrit. As Amt is decreased, the intermediate bump blends into the broad focus-point peak. 
Once again, this can be understood from Fig. 4, where we see that as Am^ is decreased, we 
require a large sensitivity to rrit in order to get an increase in C This happens at higher mo, 
and thus explains the shift in the intermediate bump to higher mo as Amt decreases. When 
the uncertainties in rrit and rrib are set to 0, we obtain a narrow peak in the focus-point 
region. This is suppressed relative to the coannihilation peak, due to the effect of the g'^ — 2 
contribution to the likelihood. 

We can now understand better Tables 1 and 2 for tan /3 = 10. For the cases with Am* 7^ 
in Table 1 and Am^ = 5 GeV in Table 2, the coannihilation peak is much higher than the 
focus-point peak, so that the 68% CL (or even the 80% CL) does not include the focus point. 
To reach the 90% CL, we need to include some part of the focus point, and this explains 
why the 68% CL is much higher than the 90% CL. The Am* = 1 GeV case in Table 2 is a 
peculiar one in which the integral over the coannihilation peak is already around 68% of the 
total integral and, because the focus point peak is flat and broad, we do not need to change 
the level much to get the 90% CL. In the cases with Am* = 0, and also Amt = 0.5 GeV 
in Table 2, the focus-point peak is also relatively high and already contributes at the 68% 
CL. Therefore we do not see an order of magnitude change between the 68% CL and the 
90% CL. 

As one would expect, the effect of the g^ — ^ constraint is more pronounced when < 0. 

This is seen in Fig. 5 for the cut with mi/2 = 300 GeV. The most starthng feature is the 

absence of the coannihilation peak at low mo when the — 2 constraint is applied. In this 

case, the focus-point region survives, because the sparticle masses there are large enough for 

the supersymmetric contribution to — 2 to be acceptably small. The broad plateau at 

intermediate mo is suppressed in this case, and the likelihood does not reach the 95% CL 

when Amt — 5 GeV. Another effect of the Higgs mass hkelihood can be seen by comparing 

the coannihilation regions for the two signs of /i when mi/2 = 300 GeV and the gfx — 2 

constraint is not applied. Because the Higgs mass constraint is stronger when // < ^, the 

^For the same Higgs mass m/i, one needs to go to a higher value of mi/2 when /i < 0. For the ehoices 
tan/3 = 10 and = 300 GeV, we find using FeynHiggs nominal values rrih = 114.1 GeV for > and 
mn = 112.8 GeV when < 0. 
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Figure 5: ^45 in Fig. 3, but for fi < and mi/2 = 300 GeV, including (excluding) the Qfj. — 2 
contribution to the global likelihood in the left (right) panel. 

coannihilation peak is suppressed when /i < relative to its height when > 0. We note 
that part of the suppression here is due to the b —>■ sj constraint, which also favours positive 

We show in Fig. 6 the likelihood function along cuts in the (mi/2,'"^o) plane for tan/3 = 
35,^0 = 0, /i < and mi/2 = 1000 GeV (left panels) and 1500 GeV (right panels). The 
Qfj, — 2 contribution to the likelihood is included in the bottom panels, but not in the top 
panels. The line styles are the same as in Fig. 3, and we note that the behaviours in the 
focus-point regions are qualitatively similar. However, at mo ~ 1000 GeV the likelihood 
function exhibits double-peak structures reflecting the locations of the coannihilation strip 
and the rapid- annihilation funnels, whose widths depend on the assumed error in mt, as can 
be seen by comparing the different line styles. 

Fig. 7 displays the likelihood function along cuts in the {mi/2,mo) plane for tan/? = 
50, Ao = 0,fi>0 and mi/2 = 800 GeV (left panels) or 1600 GeV (right panels). The 0^,-2 
contribution to the likelihood is included in the bottom panels, but not in the top panels. 
The line styles are the same as in Fig. 3, and we note that the coannihilation and focus-point 
regions even link up somewhat below the 95% CL in the case of mo = 800 GeV, if the present 
error in rrit is assumed, but only if the Qfi — 2 contribution to the likelihood is discarded. In 
this case, we can not resolve the difference between the coannihilation and funnel peaks. 
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Figure 6: As in Fig. 3 for tan/5 = 35, = 0, < and mi/2 = 1000, 1500 GeV in the 
and right panels. The — 2 constraint is included (excluded) in the bottom (top) panels 
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4 Likelihood Contours in the (mi/2,mo) Planes 



Using the fully normalized likelihood function Ctot obtained by combining both signs of 
// for each value of tan/?, we now determine the regions in the (mi/2,?Tio) planes which 
correspond to specific CLs. For a given CL, x, an iso-likelihood contour is determined such 
that the integrated volume of Ctot within that contour is equal to x, when the total volume 
is normalized to unity. The values of the likelihood corresponding to the displayed contours 
are tabulated in Table 1 (with — 2) and Table 2 (without (jf^ — 2). 

Fig. 8 extends the previous analysis to the entire (mi/2,mo) plane for tan/? = 10 and 
Aq = 0, including both signs of fi. The darkest (blue), intermediate (red) and lightest (green) 
shaded regions are, respectively, those where the likelihood is above 68%, above 90%, and 
above 95%. Overall, the likelihood for < is less than that for /j, > 0, even without 
including any information about 0^ — 2 due to the Higgs and b ^ s'j constraints. Only the 
bulk and coannihilation-tail regions appear above the 68% level, but the focus-point region 
appears above the 90% level, and so cannot be excluded. 

The highly non-Gaussian behaviour of the likelihood shown in Fig. 8 can be understood 
when comparing this figure to Fig. 3(a,b). At fixed mi/2 and for a given CL, portions of 
the likelihood function above the horizontal lines in 3(a,b) correspond to shaded regions in 
Fig. 8. The broad low-lying bump or plateau in the likelihood function at intermediate values 
of mo is now reflected in the extended features seen in Fig. 8. The extent of this plateau is 
somewhat diminished for < 0. 

The bulk region is more apparent in the left panel of Fig. 8 for > than it would 
be if the experimental error in rrit and the theoretical error in rrih were neglected. Fig. 9 
complements the previous figures by showing the likelihood functions as they would appear if 
there were no uncertainty in rrit, keeping the other inputs the same and using no information 
about Qfj, — 2. We see that, in this case, both the coannihilation and focus-point strips rise 
above the 68% CL. 

Fig. 10 is also for tan /? = 10 and Aq = 0, including both signs of /j,. This time, we include 
also the — 2 likelihood, calculated on the basis of the e^e~ annihilation estimate of the 
Standard Model contribution. This figure represents an extension of Fig. 3(c,d). In this 
case, very low values of mi/2 and mo are disfavoured when fi > 0. Furthermore, for fi < 
the likelihood is suppressed and no part of the coannihilation tail is above the 68% CL. In 
addition, none of the focus point region lies above the 90% CL for either positive or negative 
fi. However, neither of these can be excluded completely, since there are < zones within 
the 90% likelihood contour, and focus-point zones within the 95% likelihood contour. 
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Figure 8: Contours of the likelihood at the 68%, 90% and 95% levels for tan /5 = 10, = 
and fi > (left panel) or fi < (right panel), calculated using information of ruh, b S'j 
and VLcDMh^ o,nd the current uncertainties in rrit and mi,, but without using any information 
about — 




Figure 9: As in Fig. 8 hut assuming zero uncertainty in mt- 
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Figure 10: Likelihood contours as in Fig. 8, hut also including the information obtained by 
comparing the experimental measurement of g^^ — 2 with the Standard Model estimate based 
on e+e" data. 

It is important to note that the results presented thus far are somewhat dependent on the 
range chosen for mi/2; which has so far been restricted for tan /5 = 10 to < 1 TeV. In Fig. 11, 
we show the the tan/3 = 10 plane for > and yU < allowing mi/2 up to 2 TeV, including 
the gfj, — 2 constraint. Comparing this figure with Fig. 10, we see that a considerable portion 
of the focus-point region is now above the 90% CL due to the enhanced volume contribution 
at large mi/2- 

Fig. 12 is for tan/3 = 35, Aq = for both /i > and /i < 0. Fig. 13 includes the g^ — 2 
likelihood, calculated on the basis of the e~^e~ annihilation estimate of the Standard Model 
contribution, which is not included in the previous figure. In this case, regions at small mi/2 
and mo are disfavoured by the b s'-f constraint, as seen in both figures with /i < 0. At 
larger mi/2, the coannihilation region is broadened by a merger with the rapid-annihilation 
funnel that appears for large tan (3. The optional g^ — 2 constraint would prefer /x > 0, and 
in the /i < half-plane it favours larger mi/2 and mo, as seen when the left and right panels 
are compared. 

Fig. 14 is for tan /? = 50, Ao = and /i > 0. Again, the right panel includes the g^ — 2 
likelihood, calculated on the basis of the e~^e~ annihilation estimate of the Standard Model 
contribution, which is not included in the left panel. In this case, the disfavouring of regions 
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Figure 11: Likelihood contours as in Fig. 10, hut extending the range for mi 12 up to 2 TeV. 




500 1000 1500 2000 500 1000 1500 2000 

mi/2(GeV) mi/2(GeV) 

Figure 12: Likelihood contours as in Fig. 8, hut for tan/3 = 35, = and fi > (fi < 0) 
in the left (right) panel, calculated using information of rrih, b ^ s'-f and VLcDNih'^ o-nd the 
current uncertainties in rrit and mb, without the indicative information from — 2. 
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Figure 13: Likelihood contours as in Fig. 12, including the indicative information from g^j^ — 2. 

at small mi/2 and mo by the b ^ s'y constraint is less severe than in the case of tan/? = 35 
and /i < 0, but is still visible in both panels. The coannihilation region is again broadened 
by a merger with the rapid-annihilation funnel. In the absence of the (7^ — 2 constraint, both 
the coannihilation and the focus-point regions feature strips allowed at the 68% CL, and 
these are linked by a bridge at the 95% CL. However, when the optional g^ — 2 constraint is 
applied, this bridge disappears, the 90% and 95% CL strips in the focus-point region becomes 
much thinner, and the 68% strip disappears in this region. 

5 Summary 

We have presented in this paper a new global likelihood analysis of the CMSSM, incorpo- 
rating the likelihoods contributed by ruh, b ^ S'j, QcDAih'^ and (optionally) g^ — 2. We have 
discussed extensively the impacts of the current experimental uncertainties in mt and mj,, 
which affect each of rrih, b ^ s'j and QcDAih'^- In particular, the widths of the coannihilation 
and focus-point strips are sensitive to the uncertainties in mt and m^, and a low-lying plateau 
in the likelihood is found with the present uncertainty Amt = 5 GeV. 

We recall that the absolute values of the likelihood integrals are not in themselves mean- 
ingful, but their relative values do carry some information. Generally speaking, the global 
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Figure 14: Likelihood contours as in Fig. 12, hut for tan/? = 50, = and > 0, without 
(left panel) and with (right panel) the indicative information from — 2. 

likelihood function reaches higher values in the coannihilation region than in the focus-point 
region, as can be seen by comparing the entries with and without parentheses in Table 3. 
This tendency would have been reversed if the uncertainty in rrit had been neglected, as seen 
in Table 4, but the preference for the coannihilation region is in any case not conclusive. 

Table 3 also displays the integrated likelihood function for different values of tan jS and 
the sign of /x, exhibiting a weak general preference for /i > if the g^ — 2 information is used. 
If this information is not used, /i < is preferred for tan/? = 35, but ;U > is still preferred 
for tan (3 = 10. There is no significant preference for any value between tan (3 = 10 and the 
upper limits ^ 35 and ^ 50 where electroweak symmetry breaking ceases to be possible in 
the CMSSM, though we do find a weak preference for tan/? = 50 and /i > 0. 

In the foreseeable future, the analysis in this paper could be refined with the aid of 
improved measurements of rrit at the Fermilab Tevatron collider, by refined estimates of 
nih, by better determinations of ^cDuh'^ and more experimental and theoretical insight into 
g^ — 2, in particular. One could also consider supplementing our phenomenological analysis 
with arguments based on naturalness or fine-tuning, which would tend to disfavour larger 
values of mi/2 and mo. However, in the absence of such theoretical arguments, our analysis 
shows that long strips in the coannihilation and focus-point regions cannot be excluded on 
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Table 3: Integrals of the global likelihood function integrated over the (mi/2,mo) planes for 
various values of tan f3, in the coannihilation and funnel (focus-point) region. We use the 
range mi/2 < 'iTeV , except for the second row of tan (5 = 10 case, where the range mi/2 < 
1 TeV is used. 





incl. g^-2 


without 


^?M-2 


tan P 


H> 


jj, <0 


/X > 




10 


41.7 (5.9) 


2.1 (4.8) 


2329 (1052) 


1147 (984) 




11.7 (2.9) 


2.1 (1.7) 


2329 (176) 


1117 (387) 


35 


33.9 (11.6) 


25.9 (5.5) 


1428 (1596) 


8690 (1270) 


50 


231.9 (6.84) 




13096 (1117) 





Table 4: As in Table 3, but assuming zero uncertainty in rut- 



tan P 


incl. g^-2 


without g^ — 2 




/X < 


fi> 


/i < 


10 


44.9 (69.1) 
44.9 (21.4) 


2.6 (67.7) 
2.6 (20.2) 


2425 (12916) 
2425 (3922) 


1485 (13442) 
1485 (4144) 


35 


33.5 (90.9) 


26.4 (58.3) 


1451 (15377) 


8837 (12589) 


50 


195.0 (60.4) 




13877 (10188) 
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the basis of present data. The preparations for searches for supersymmetry at future coUiders 
should therefore not be restricted to low values of mi/2 and mo- 
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